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FOREWORD 


The present report is one of a series of six reports, published simul- 
taneously, which describe analyses and computational procedures for: 1) pre- 

diction of the in-depth response of charring ablation materials, based on one- 
dimensional thermal streamtubes of arbitrary cross-section and considering 
general surface chemical and energy balances, and 2) nonsimilar solution of 
chemically reacting laminar boundary layers, with an approximate formulation 
for unequal diffusion and thermal diffusion coefficients for all species and 
with a general approach to the thermochemical solution of mixed equilibrium- 
nonequilibrium, homogeneous or heterogeneous systems. Part I serves as a 
summary report and describes a procedure for coupling the charring ablator 
and boundary layer routines. The charring ablator procedure is described in 
Part II, whereas the fluid-mechanical aspects of the boundary layer and the 
boundary-layer solution procedure are treated in Part III. The approximations 
for multicomponent transport properties and the chemical state models are 
described in Parts IV and V, respectively. Finally, in Part VI an analysis • 
is presented for the in-depth response of charring materials taking into ac- 
count char-density buildup near the surface due to coking reactions in depth. 

The titles in the series are: 


Part I Summary Report: An Analysis of the Coupled Chemically Reacting 

Boundary Layer and Charring Ablator, by R. M. Kendall, E. P. 
Bartlett, R. A. Rindal, and C. B. Moyer. 

Part II Finite Difference Solution for the In-depth Response of Charring 
Materials Considering Surface Chemical and Energy Balances, by 
C. B. Moyer and R. A. Rindal. 


Part III Nonsimilar Solution of the Multicomponent Laminar Boundary Layer 

by an Integral Matrix Method, by E. P. Bartlett and R. M. Kendall. 

Part IV A Unified Approximation for Mixture Transport Properties for Multi- 
component Boundary-Layer Applications, by E. P. Bartlett, R. M. 
Kendall, and R. A. Rindal. 


Part V A General Approach to the Thermochemical Solution of Mixed Equilib- 
rium-Nonequilibrium, Homogeneous or Heterogeneous Systems, by 
R. M. Kendall. 


Part VI An Approach for Characterizing Charring Ablator Response with In- 
depth Coking Reactions, by R. A. Rindal. 


This effort was conducted for the Structures and Mechanics Division of 
the Manned Spacecraft Center, National Aeronautics and Space Administration 
under Contract No. NAS9-4599 to Vidya Division of Itek Corporation with Mr. 
Donald M. Curry and Mr. George Strouhal as the NASA Technical Monitors. The 
work was initiated by the present authors while at Vidya and was completed 
by Aerotherm Corporation under subcontract to Vidya (P.0. 8471 V9002) after 
Aerotherm purchased the physical assets of the Vidya Thermodynamics Depart- 
ment. Dr. Robert M. Kendall of Aerotherm was the Program Manager and Prin- 
cipal Investigator. 
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ABSTRACT 

Self-consistent approximations for binary diffusion coefficients, multi- 
component thermal diffusion coefficients and mixture viscosity and thermal 
conductivity are presented which greatly simplify the solution of multicom- 
ponent laminar boundary-layer problems. The basic premise of the approach is 
a bifurcation of binary diffusion coefficients, 2K ^ « D/F^F^ , where D is 
a reference diffusion coefficient and F^ is a diffusion factor for species 
i . This expression is exact for binary or ternary systems but is approximate 
for larger systems since v(v - l)/2 are replaced by v F^ where v is 

the number of species. The F^ are determined by a least-squares correlation 
of actual diffusion coefficient data. The adequacy of the correlation is 
demonstrated for a number of chemical systems. 

The primary utility of the bifurcation approximation in multicomponent 
boundary-layer applications is that it enables explicit solution of the 
Stefan -Maxwell relations for the diffusive mass fluxes, j ^ . Use of these 
relations, in turn, permits the Shvab-Zeldovich transformation of the species 
conservation equations without introduction of the concentration-dependent 
multicomponent diffusion coefficients, . In addition, an approximation 

to the modified Stefan -Maxwell relations provides the framework for the 
rationalization of a film coefficient model for diffusional heat and mass 
transfer . 

The approximation for multicomponent thermal diffusion coefficients is 
based upon generalization of a semi-empirical equation for correlating binary 
thermal diffusion data. The basic hypothesis in this generalization to multi- 
component systems is that the thermal diffusion of each species behaves as if 
it were in a binary mixture of that species and a species representative of 
the mixture as a whole. The result is algebraically very simple and enables 
consideration of multicomponent thermal diffusion and diffusion thermo in large 
complex systems with negligible additional effort when done in conjunction with 
the approximate formulation for unequal binary diffusion coefficients. An 
evaluation of the accuracy of the thermal diffusion approximation in multi- 
component systems has not been attempted to date. 

Finally, approximate relations for mixture viscosity and thermal conduc- 
tivity of the Suther land-wassili jewa type are simplified by use of the 
bifurcation approximation for binary diffusion coefficients. 
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A UNIFIED APPROXIMATION FOR MIXTURE TRANSPORT 
PROPERTIES FOR MULTICOMPONENT BOUNDARY- LAYER APPLICATIONS 


SECTION 1 
INTRODUCTION 


Consideration of unequal diffusion and thermal diffusion coefficients 
substantially complicates the solution of multicomponent laminar boundary- 
layer problems when these transport properties are based on rigorous kinetic 
theory. For this reason, diffusion coefficients have often been assumed equal 
and thermal diffusion has been neglected in many situations where these approx- 
imations can in fact produce substantial errors, a correlation for binary 
diffusion coefficients, first considered by Bird,'*' is utilized herein which 
permits consideration of unequal diffusion coefficients for all species in a 
relatively precise manner with little increase in algebraic complexity and 
computational time relative to assumed equal diffusion coefficients. This 
correlation also provides the framework for an approximation for multicomponent 
thermal diffusion coefficients which permits thermal diffusion to be conven- 
iently taken into consideration. The correlation for binary diffusion coef- 

2 

ficients is also employed to simplify Sutherland-Wassil jewa -type approximate 
expressions for mixture viscosity and thermal conductivity, thus completing 
a unified theory for multicomponent boundary-layer transport properties. The 
correlation for binary diffusion coefficients is described in Section 2. The 
approximate relations for mixture viscosity, mixture thermal conductivity, 
and multicomponent thermal diffusion coefficients are developed in Sections 3 
through 5, respectively. 


SECTION 2 

A BIFURCATION APPROXIMATION FOR BINARY 
DIFFUSION COEFFICIENTS 

There have been few attempts to solve the multicomponent boundary-layer 
equations for chemical systems involving large numbers of chemical species 
while considering unequal diffusion coefficients for all species. The primary 
hindering factor has been the complex interdependence of the mass diffusion 
flux, j^, upon the species concentrations and concentration gradients. An 
approximation for binary diffusion coefficients is introduced herein which 
permits substantial simplification of the expression for j ^ with attendant 
simplification of multicomponent boundary- layer problems. 


In the past the conventional approach has been to assume equal diffusion 
coefficients for all species. When this is done, the diffusion mass flux of 
species i, j^, is related explicitly to its mass fraction gradient by Fick*s 
law 

dK. 

- - »-»i2 nr w 


where t * ie kinary diffusion coefficient for the system, p is the 

density, and is the mass fraction of species i. This has merit for some 

problems when it can be argued, for example, that all light species are dif- 
fusing in the same direction. However, in general it is a poor approximation 
and is justifiable only as a convenience. When unequal diffusion coefficients 

have been taken into consideration, the number of chemical species has generally 

3 

been limited to a relatively small number. For example, Scala and Gilbert 
considered a nine-component model for the problem of graphite ablation in air. 

A chemical system of this size does not allow consideration of trace species 
(these can be important, for example, in radiative heating and reentry com- 
munication problems) . Furthermore, it is generally inadequate for more com- 
plicated materials such as reinforced organics. 

In the absence of diffusion due to pressure gradients and external forces, 

2 

the mass diffusion fluxes are given implicitly by the Stefan-Maxwell relations. 
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where x^ is the mole fraction of species i, T is the temperature, iK ^ is 
the binary diffusion coefficient for species i and j, and is the multi- 

component thermal diffusion coefficient for species i. Utilization of these 
relations in conjunction with the boundary- layer conservation equations is 
awkward even in the absence of thermal diffusion effects as a result of the 
implicit behavior of j ^ on mole fractions and their gradients. An explicit 
relation for the mass diffusion flux can be obtained by employing multicompo- 
nent diffusion coefficients, D. . (see Ref. 2) : 
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where ! Yf { . and T7\ are the molecular weights of species i and of the system. 
Unfortunately, each of the (v -1) v depend upon local concentrations and 

upon (v - 1) v/2 symmetric binary diffusion coefficients, where v is the 
total number of species being considered. Equation (3) has generally been 
used in preference to Eq. (2) for multicomponent boundary- layer problems 
(e.g.. Refs. 3, 4 and 5). 

The reason that the Stefan-Maxwell relations (Eq. (2)) cannot be arranged 
into an explicit relationship for j ^ without introduction of the is 

that the contributions of species i and j to the binary diffusion coefficient 

are inseparable. According to first-order kinetic theory, binary diffusion 
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coefficients are expressible as 
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where 


T temperature, °K 

77(^,771 j molecular weights of species i, j 
P pressure in atmospheres 

ttg . . collision cross section £or "rigid sphere" molecule i 
ID o 

interacting with j , A 


^iV^* integral expression which depends upon the particular inter- 
molecular potential function which is utilized. 

In Ref. 1, Bird showed that a bifurcation of the effects of species i 
and j in Eq. (4) permits explicit solution of the Stefan-Maxwell relations 
for j ^ in terms of gradients and properties of species i and of the system 
as a whole. This approximation is utilized herein in the form 





(5) 


where D is a reference diffusion coefficient and F^ might be termed a 
diffusion factor for species i. Equation (5) is exact for a ternary system 
but approximate for larger systems, since (v - 1) v/2 diffusion coefficients, 
j, are replaced by v diffusion factors, F^. Equation (5) should be 
looked upon as a correlation equation for actual binary diffusion coefficient 
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data. This correlation can be obtained by various means. The present ap- 
proach has been to first define 
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2.6280 x 10" 3 
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with T in °K, P in atmospheres, and a in A. (The subscript "ref ,r refers 
to a reference species (often but conceivably fictional) . D is thus the 
self-diffusion coefficient of that species.) The are then determined by 

a least-squares correlation of the for all diffusing pairs in the chemi- 

cal system of interest. With this approach the pressure dependence and the 
majority of the temperature dependence of are absorbed into the D so 

that the are independent of pressure and have only a secondary tempera- 

ture dependence. All species have been treated on an equal basis (i.e., 
without the use of any weighting factors) . This has the decided advantage 
that the F^ are then independent of concentration and can be determined 
a priori for a given set of chemical species. The accuracies which have been 
obtained using this procedure have been surprisingly good as will be shown 
later. 


The development of the explicit relation for j ^ from the Stefan-Maxwell 
relations (Eq. (2)) is presented in Appendix A. The result can be expressed 
as 
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The Z i is a quantity for species i which sums to unity and lies between a 
mass fraction and a mole fraction, whereas |i , \i and are properties 

of the system as a whole. It will be shown that it may often be consistent 
with the level of the approximation to consider the independent of tem- 

perature. In this case Eq. (7) reduces to 
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j + fi-12 

-’i T 3y 


pDta 2 
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K i } ^2 
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( 12 ) 


In the absence of thermal diffusion and assuming dz^/dy to be large 
compared to £(z^ - K^) /m 2 U 9 Eq- (12) simplifies to 


pdh 2 az ± 


(13) 


The accuracy of this approximation is assessed in Appendix B. The approxi- 
mate Eq. (13) provides the basis for rationalizing a film-coefficient rela- 
tion for diffusional mass transfer in the multicomponent boundary layer. 

This model is described in Appendix B and represents an extension of the 
model presented in Ref. 7. 

The choice of Eq. (11) , (12) , or (13) depends upon whether or not the 

additional approximations embodied in the latter two equations are satisfac- 
tory for the particular chemical system under consideration and the desired 
degree of accuracy. Equation (12) is recommended for boundary- layer calcu- 
lations and is utilized in Ref. 6. Equation (13) is appropriate for heat 
and mass transfer calculations utilizing overall convective transfer coeffi- 
cients . 

The accuracy of the correlation (Eq. (5)) was investigated by Bird 1 for 
a five-component mixture of N 2 , 0 2# C0 2# CO, and H 2 and shown to be surpris- 
ingly good, the maximum error in any ih ^ being 4 percent. In order to 
establish more generally the adequacy of the approximation, correlations 
have been performed for several larger chemical systems. In order to inter- 
pret the results of these correlations, it is pertinent to first briefly 
review elementary procedures for predicting binary diffusion coefficients. 
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The kinetic theory relation for was presented as Eq. (4). In the 

E rigid sphe 

is unity by definition. Since T and P are thermodynamic proper- 


simplest and crudest of approximations, the model of rigid spheres, the 

of 1 * 1 )* 

ties and the V[. are known constants, the only additional information required 


is 

law 
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For the rigid sphere model, this is given by an additive combining 
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and 


are determined empirically or estimated for each molecu- 


where a ^ and 
lar species i and j individually. 

There are several improved models for the intermolecular potential func- 
tion with additional empirical constants. Probably the most popular is the 

(i ,*)* 


for 


Lennard-Jones model which has two empirical constants. The ft: 

2 

this model are tabulated in terms of a reduced temperature, T^ , defined as 


T h 


T 

e ij A 


(15) 


is the maximum energy of attraction of the two molecules i and 

The e^/k is the order of 100 for 


where 

j, and k is the Boltzmann constant. The e . ./k 

1 \ 

most simple molecules and ft^' * varies from 1.44 for T* = 1.0 to 0.42 
for T* = 400. The o . A and e,._./k are not well known for many interactions 

a .. ^ and 


ID i r 

so they are often estimated by Eq. (14) for 


iD 


'ID 
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(16) 


The Lennard-Jones model is not always accurate, especially for nonsymmetric 
molecules (such as polarized molecules or long chains) and at high tempera- 
tures (where ionization and excitation occur) . On the other hand, this model 
is useful for preliminary calculations of complex chemical systems, since the 
two empirical constants can, in the absence of experimental data, be estimated 
by various correlation techniques, a situation not always possible with more 
sophisticated models. 

Results of a correlation for a 9-component (36 O-N-C system, rep- 

resentative of the boundary layer over graphite ablating into air, are pre- 
sented in Table I. The F^ were obtained by employing a least-squares fit 
to j for all diffusing pairs based on the Lennard-Jones potential using 
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TABLE I.- CORRELATION OF BINARY DIFFUSION COEFFICIENTS FOR AN 
OXYGEN-NITROGEN-CARBON SYSTEM BASED ON THE LENNARD-JONES 
POTENTIAL WITH FORCE DATA FROM SVEHLA (Ref. 8) . 

TEMPERATURE = 12,000°R, PRESSURE - 1 ATM. 


SPECIES 

FROM 

p i to 

FROM PRES- 

ERROR USING 

ERROR IF ALL 



KINETIC THEORY 

CORRELAT I ON 

PRESENT 

ARE AS- 
SUMED EQUAL 



2 

0 


CORRELATION 

i 

j 

(ft /sec) x 100 

(ft /sec) x 100 

(PERCENT) 

(PERCENT) 

0 

02 

5.6458 

0.7399' 

5.5575 

- 1.6 

- 16.8 

0 

N 

7.3372 


7.5274 

2.6 

- 36.0 

0 

N2 

5.3995 


5.3837 

- 0.3 

- 13.0 

0 

CO 

5.4662 


5.4379 

- 0.5 

- 14.1 

0 

C02 

4*4638 


4.3762 

- 2.0 

5.2 

0 

C 

8.0754 


8.3663 

3.6 

- 41.9 

0 

C3 

5.1820 


5.0859 

- 1.9 

- 9.4 

0 

CN 

5.3620 


5.3697 

0.1 

- 12.4 

02 

N 

5.6566 

1.0000 

5.5695 

- 1.5 

- 16.9 

02 

N2 

3.9611 


3.9834 

0.6 

18.6 

02 

CO 

4.0028 


4.0235 

0.5 

17.3 

02 

C02 

3.1637 


3.2380 

2.3 

48.4 

02 

C 

6.3129 


6.1902 

- 1.9 

- 25.6 

02 

C3 

3.7100 


3.7630 

1.4 

26.6 

02 

CN 

3.9623 


3.9731 

0.3 

18.4 

N 

N2 

5.4277 

0.7383 

5.3953 

- 0.6 

- 13.5 

N 

CO 

5*^763 


5.4496 

- 0.5 

- 14.2 

N 

C02 

4.5136 


4.3857 

- 2.8 

4.0 

N 

C 

7.9727 


8.3844 

5.2 

- 4 1 . 0 

N 

C3 

5.2069 


5.0969 

- 2.1 

- 9.8 

N 

CN 

5.3784 


5.3813 

0. 1 

- 12.7 

N2 

CO 

3.8943 

1.0323 

3.8977 

0.1 

20.6 

N2 

C02 

3.1114 


3.1367 

0.8 

51.0 

N2 

C 

6.0528 


5.9967 

- 0.9 

- 22.4 

N2 

C3 

3.6214 


3.6454 

0.7 

29.7 

N 2 

CN 

3.8603 


3.8488 

- 0.3 

21.7 

CO 

C02 

3.1390 

1.0220 

3.1683 

0.9 

49.7 

CO 

C 

6.II84 


6.0570 

- 1.0 

- 23.2 

CO 

C3 

3.6584 


3.6820 

0.6 

28.4 

CO 

CN 

3.8938 


3.8875 

- 0.2 

20.6 

C02 

C 

4.9902 

1.2700 

4.8745 

- 2.3 

- 5.9 

C02 

C3 

2.8753 


2.9632 

3.1 

63*4 

C02 

CN 

3.1245 


3.1286 

0.1 

50.3 

C 

C3 

5.7767 

0.6643 

5.6649 

* 1.9 

- 18.7 

C 

CN 

6.0033 


5.9811 

- 0.4 

- 21.8 

C3 

CN 

3.6276 

1.0927 

3.6359 

0.2 

29.5 

CN 



1.0350 




AVERAGE 

ABSOLUTE ERROR 



1.3 

24.2 
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molecular properties suggested in Ref. 8. Utilizing these F^, diffusion 
coefficients were computed from Eq. (5) and compared to the original j . 

It can be seen that substantial improvement over the equal diffusion coeffi- 
cient model is obtained, the average absolute error in 3 ^ being reduced 
by more than an order of magnitude. Even more significant, however, is the 
comparison with the original 2K average absolute error being 1.3 per- 
cent and the maximum error in any one single 3 ^ 5.2 percent. Better accu- 

racies than this are probably of academic interest in view of the uncertain- 
ties in the diffusion coefficients themselves. 

Although this correlation is very encouraging, it was obtained with 
consistent use of a single, admittedly crude model and thus may not be a 
severe test of the correlation procedure. Therefore, a 6-component (15 JK j) 
hydrogen-oxygen system was investigated for two different sets of 3 * . : (1) data 

1 J Q 

based on the Lennard-Jones model and (2) data recently compiled by Svehla 
based on improved models for the intermolecular potentials for the H-O system. 
The factors and the resulting errors in are presented for these 

two sets of diffusion data in Table II. The more reliable set of diffusion 
data (Table II (b) ) is seen to correlate better than do the data where the 
Lennard-Jones potential and combining rules (Eqs. (14) and (16)) are employed 
consistently (Table 11(a)), the average absolute error being reduced from 10.8 
to 4.8 percent and the maximum error in any one diffusion coefficient being 
reduced from 47.6 to 11.3 percent. The primary reason for this is that the 
improved 3 ^ for H - (based upon the recent high- temperature calcula- 
tions of Ref. 10) is higher and thus more in line with the value required for 
a good least-squares fit. The correlation of Table 11(b) is considered reason- 
ably good in view of the molecular weight variation of 32 to 1. 

Results for a 16-component (120 ^j) 0-N-C-H system, representative of 
the boundary layer over graphite-phenolic ablating in air, are presented in 

9 

Table III. These were based on the procedures suggested by Svehla. 

In this system, 98 of the 120 are represented within 5 percent, the lar- 

gest single error is 16.7 percent and the average error is 3.7 percent. Again 
the correlation is deemed successful considering that 120 numbers are repre- 
sented by 16 numbers and that the molecular weight variation is 44 to 1. 

One set of diffusion coefficient data was considered which did not cor- 
relate well, namely, the data utilized in Ref. 3 for the same 0-N-C system 
studied in Table I. The results of the correlation are presented in Table IV. 
In this case there is about a 25-percent improvement over assumed equal dif- 
fusion coefficients rather than the usual order-of-magnitude improvement. 

The explanation for the lack of correlation is as follows: the Lennard-Jones 

model was used for interactions involving two polyatomic species, but the 
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TABLE II.- CORRELATION OF BINARY DIFFUSION COEFFICIENTS FOR A 
HYDROGEN-OXYGEN SYSTEM USING PRESENT METHOD. 

TEMPERATURE - 12,000°R, PRESSURE - 1 ATM. 

(a) Diffusion coefficients calculated using Lennard-Jones 
potential with force data from Svehla (Ref. 8) 


species 

FROM 

p i 

i>i- 

; FROM PRES- 

ERROR USING 

ERROR IF AL 



KINETIC THEORY 


ENT 

CORREL AT I ON 

PRESENT 

ARE AS 



0 



0 

CORRELATION 

SUM ED F 0 U A 

i 

j 

(ft /sec) x 100 

(ft /sec) x 100 

(PERCENT) 

( PERCENT ) 

H 

H2 

36.0260 

0.24713 


53. 1613 

47.6 

- 63. 1 

H 

H2U 

23.9891 



23.7639 

- 8.6 

- 48.9 

H 

0 

26.6238 



24.7360 

- 7.1 

- 50.1 

H 

02 

22.8038 



19.7757 

-13.3 

- 41.7 

H 

OH 

26.4341 



24.3147 

- 8.0 

- 49.7 

H2 

H20 

17.3862 

0.3720 


15.7877 

- 9.2 

- 23.5 

H2 

0 

17.7166 



16.4335 

- 7.2 

- 24.9 

H2 

02 

15.0083 



13.1381 

-12.5 

- 11.4 

H2 

OH 

17.5759 



16.1537 

- 8. 1 

- 24.4 

H2U 

0 

7.0928 

0.8322 


7.3461 

3.6 

87.4 

H20 

02 

5.2795 



5.8730 

11.2 

151.6 

H20 

OH 

6.9078 



7.2210 

4.5 

92.4 

0 

02 

5.6458 

0.7995 


6.1132 

8.3 

135.4 

CJ 

OH 

7.2643 



7.5163 

3.5 

82.9 

02 

OH 

5.4946 

1.0000 


6.0091 

9.4 

141.9 

OH 



0.8133 





AVERAGE 

ABSOLUTE ERROR 




10.8 

68.6 


(b) Diffusion coefficients calculated using values for collision 
cross-sections suggested by Svehla (Ref. 9) 


H 

H2 

67.6000 

0.2208 

74.4024 

10.1 

- 77.1 

H 

H20 

28.3200 


27.0030 

- 4.7 

^ 45.4 

H 

0 

27.7200 


30.8482 

11.3 

- 44.3 

H 

02 

24.5500 


22.5734 

-8.1 

- 37.1 

H 

OH 

29.5900 


27.5549 

- 6.9 

- 47.8 

M2 

H20 

19.5800 

0.3034 

19.6568 

0.4 

- 21.1 

H2 

0 

23.6000 


22.4560 

- 4 • H 

- 34.5 

H2 

02 

17. 1900 


16.4323 

- 4.4 

- 10.1 

H2 

OH 

20. 1600 


20.0586 

- 0.5 

- 23.3 

H2G 

0 

8.2950 

0.8360 

8. 1500 

- 1.7 

86. 3 

H20 

02 

5.7150 


5.9638 

4.4 

170.4 

H2U 

OH 

7.1450 


7.2799 

1.9 

116.3 

0 

02 

6.8500 

0.7317 

6.8131 

- 0.5 

125.6 

0 

OH 

8.6060 


8.3166 

- 3.4 

79.6 

02 

OH 

5.5520 

1.0000 

6.0857 

9.6 

178.3 

(JH 



0.8192 




AVERAGE 

ABSOLUTE ERROR 



4.8 

73. 1 
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TABLE III.- CORRELATION OF BINARY DIFFUSION COEFFICIENTS FOR AN 
OXYGEN-NITROGEN-CARBON-HYDROGEN SYSTEM BASED ON DATA OF SVEHLA 
(Refs. 8 and 9). TEMPERATURE » 12,000°R, PRESSURE - 1 ATM. 


species 

•^ij FROM Fi 

•^ij FROM PRES- 

ERROR USING 

ERROR IE ALL 



KINETIC THEORY 

ENT CORRELATION 

PRESENT 

«^ij ARE AS- 



0 

2 . 

CORRELATION 

SUMED EOUAL 

i 

j 

(ft /sec) x 100 

(ft /sec) x 100 

(PERCENT) 

( PERCENT) 

0 

02 

6,8500 0,7393 

6,0528 

-11.6 

- 4.5 

0 

N 

7.3372 

7.6554 

4.3 

- 10.0 

0 

N2 

5.3995 

5.6277 

4.2 

21*1 

0 

CO 

5.4662 

5.6846 

4.0 

19.7 

0 

C02 

4.4638 

4.6277 

3.7 

46.5 

0 

C 

8.0754 

8.3865 

3.9 

- 19.0 

0 

C 3 

5.1820 

5.3597 

3.4 

26.2 

0 

CN 

5.3620 

5.5958 

4.4 

22.0 

0 

H 

27.7200 

29.8130 

7.6 

- 76.3 

0 

H2 

23.6000 

20.4311 

-13.4 

- 72.2 

0 

H20 

8.2950 

7.5057 

- 9.5 

- 21.1 

0 

OH 

8.6060 

7.7923 

- 9.5 

- 24.0 

0 

CH4 

5.8190 

6.0648 

4,6 

12.4 

0 

C2H 

4.8947 

5,0977 

4.1 

33.6 

0 

HCN 

4.8625 

5.0401 

3.7 

34.5 

02 

N 

5.6566 1.0000 

5.6595 

0.1 

15.6 

02 

N2 

3.9611 

4.1604 

5.0 

65.1 

02 

CO 

4.0028 

4.2025 

5.0 

63*4 

02 

C02 

3.1637 

3.4212 

8.1 

106*7 

02 

C 

6.3129 

6.2000 

- 1.8 

3*6 

02 

C3 

3.7100 

3.9624 

6.8 

76.3 

02 

CN 

3.9623 

4.1369 

4.4 

65.1 

02 

H 

24.5500 

22.0402 

-10.2 

- 73*5 

02 

H2 

17.1900 

15.1043 

-12.1 

- 62*0 

02 

H20 

5.7150 

5.5489 

- 2.9 

14.5 

02 

OH 

5.5520 

5.7607 

3.8 

17.8 

02 

CH4 

4.4735 

4.4984 

0.6 

46.2 

02 

C2H 

3.6310 

3.7686 

3.8 

80.1 

02 

HCN 

3.5678 

3.7260 

4.4 

03.3 

N 

N2 

5.4277 0.7907 

5.2620 

- 3.1 

20.5 

N 

CO 

5,4763 

5.3153 

- 2.9 

19.4 

N 

C02 

4.5136 

4.3270 

- 4.1 

44.9 

N 

C 

7.9727 

7.8416 

- 1*6 

- 17.9 

N 

C3 

5.2069 

5.0115 

- 3.8 

25.6 

N 

CN 

5.3784 

5.2323 

- 2.7 

21.6 

N 

H 

25.5139 

27.8700 

9.3 

- 74.4 

N 

H2 

17.1218 

19,1036 

11.6 

- 61.8 

N 

H20 

6.9743 

7.0181 

0.6 

- 6.2 

N 

OH 

7.1732 

7.2861 

1.6 

- 8.8 

N 

CH4 

5.7836 

5.6895 

- 1.6 

13*1 

N 

C2H 

4.9083 

4.7665 

- 2.9 

33*3 

N 

HCN 

4.8645 

4.7126 

- 3.1 

34.5 
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TABLE III.- Continued 


species 

FROM 

F i 

FROM PRES- 

ERROR USING 

ERROR IF ALL 



KINETIC THEORY 

ENT 

CORREL AT ION 

PRESENT 

ARE AS- 






CORRELATION 

SUMED EQUAL 

i 

j 

(ft 2 /sec) x 100 

(ft /sec) x 100 

< PERCENT ) 

(PERCENT) 

N2 

CO 

3 • 8 943 

1.0756 

3.9074 

0.3 

68.0 

N 2 

C02 

3.1114 


3.1809 

2.2 

110.2 

N2 

C 

6.0528 


5.7645 

- 4.8 

8.1 

N2 

C 3 

3.6214 


3.6840 

1.7 

80.6 

N2 

CM 

3.8603 


3.8463 

- 0.4 

69.4 

N2 

H 

21.3750 


20.4922 

- 4.1 

- 69.4 

N2 

H2 

14.1671 


14.0435 

- 0.9 

- 53.8 

N2 

H20 

5.0300 


5.1591 

2.6 

30.0 

N? 

OH 

5.2629 


5.3561 

1.8 

24.3 

N2 

CH4 

4.3182 


4.1624 

- 3.1 

51.5 

N2 

C2H 

3.5367 


3.5039 

- 0.9 

84.9 

N2 

hcn 

3.4655 


3.4643 

- 0.0 

84.8 

CO 

C02 

3.1390 

1.0647 

3.2131 

2.4 

99.2 

CO 

C 

6.1134 


5.8229 

- 4.8 

9.9 

CO 

C 3 

3.6504 


3.7213 

1.7 

72.0 

CO 

CN 

3.8938 


3.8853 

- 0.2 

64.8 

CO 

H 

21.6122 


20.6996 

- 4.2 

- 66.9 

CO 

H2 

14.2296 


14.1856 

- 0.3 

- 50.4 

CO 

H20 

5.1001 


5.2113 

2.2 

27.1 

CO 

OH 

5.3305 


5.4103 

1.5 

22.3 

CO 

CH4 

4.3595 


4.2248 

- 3.1 

51.5 

CO 

C2H 

3.5680 


3.5394 

- 0.8 

80.9 

CO 

HCN 

3.5023 


3.4994 

- 0.1 

62.9 

C02 

C 

4.9902 

1.3079 

4.7402 

- 5.0 

35.0 

C02 

C 3 

2.8753 


3.0294 

5.4 

111.3 

C02 

CN 

3.1245 


3.1629 

1.2 

102.4 

C02 

H 

18.4881 


16.8510 

- 8.9 

- 59.4 

C02 

H2 

12.2917 


11.5481 

- 6.0 

- 39.1 

C02 

H20 

4.1217 


4.2424 

2.9 

58.7 

C02 

OH 

4.3441 


4.4044 

1.4 

50.6 

CC2 

CH4 

3.5835 


3.4393 

- 4.0 

82.5 

C02 

C2H 

2.8685 


2.8613 

0.4 

128.0 

C02 

HCN 

2.7965 


2.8488 

1.9 

133.9 

C 

C 3 

5.7767 

0.7219 

5.4901 

- 5.0 

13.2 

c 

CN 

6.0033 


5.7319 

- 4.5 

8.9 

c 

H 

26.1719 


30.5380 

16.7 

- 75.0 

c 

H2 

18.0635 


20.9280 

15.9 

- 63.8 

c 

H20 

7.5630 


7.6883 

1.7 

- 13.5 

c 

OH 

7.9334 


7.9818 

0.6 

- 17.5 

c 

CH4 

6.3330 


6.2328 

- 1.6 

3.3 

c 

C2H 

5.3831 


5.2217 

- 3.0 

21.5 

c 

HCN 

5.3406 


5.1626 

- 3.3 

22.5 
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TABLE III.- Concluded 


SPECIES 

FROM 

F . 
1 

hi 

FROM PRES- 

ERROR USING 

ERROR IF ALL 



KINETIC THEORY 

ENT 

CORRELATION 

PRESENT 

ARE AS- 
SUMED EQUAL 



2 



2 

CORRELATION 

i 

j 

(ft /sec) x 100 


(ft 

/sec) x 100 

( PFRCENT ) 

< PERCENT) 

C3 

CN 

3.6276 

1.1293 


3.6632 

1.0 

80.3 

C3 

H 

21.0069 



19.5166 

- 7.1 

- 68.9 

C3 

H2 

13.9792 



13.3749 

- 4.3 

- 53.2 

C 3 

H20 

4.8271 



4.9135 

1.8 

35.5 

C3 

OH 

5.0416 



5.1011 

1.2 

29.7 

C 3 

C H4 

4.1210 



3.9833 

- 3.3 

50.7 

C3 

C2H 

3.3265 



3.3371 

0.3 

96.6 

C 3 

HCN 

3.2583 



3.2994 

1.3 

100.7 

CN 

H 

20.9403 

1.0817 


20.3763 

- 2.7 

- 68.8 

CM 

H2 

13.8853 



13.9641 

0.6 

- 52.9 

CN 

H 20 

4.9948 



5.1300 

2.7 

30.9 

CN 

OH 

5.2299 



5.3259 

1.8 

25.1 

CN 

CHA 

4.2953 



4.1506 

- 3.2 

52.3 

CN 

C2H 

3.5330 



3.4841 

- 1.4 

85.1 

CN 

HCN 

3.4626 



3.4447 

- 0.5 

88.9 

H 

H2 

67.6000 

0.2030 


74.3967 

10.1 

- 90.4 

H 

H20 

28.3200 



27.3309 

- 3.5 

- 76.9 

H 

OH 

29.5900 



28.3745 

- 4.1 

- 78.0 

H 

CH4 

20.3467 



22.1568 

8.9 

- 68.0 

H 

C2H 

18.6611 



18.5624 

- 0.5 

- 64.9 

H 

HCN 

18.8560 



18.3525 

- 2.7 

- 65.3 

H2 

H20 

19.5000 

0.2963 


18.7301 

- 4.3 

- 66.7 

H2 

OH 

20.1600 



19.4453 

- 3.5 

- 67.6 

H2 

CH4 

13.7590 



15.1843 

10.4 

- 52.4 

H2 

C2H 

12.5045 



12.7210 

1.7 

- 47.7 

H2 

HCN 

12.5953 



12.5771 

- 0.1 

- 48.1 

H20 

OH 

7.1450 

0.8064 


7.1436 

- 0.0 

- 8.5 

H20 

CH4 

5.4665 



5.5782 

2.0 

19.6 

H20 

C2H 

4.5559 



4.6733 

2.6 

43.6 

H20 

HCN 

4.5242 



4.6205 

2.1 

44.6 

OH 

CH4 

5.6987 

0.7767 


5.7912 

1.6 

14.8 

OH 

C2H 

4.7817 



4.8517 

1.5 

36.8 

OH 

HCN 

4.7308 



4.7969 

1.2 

38.0 

CH4 

C2H 

3.9244 

0.9948 


3.7806 

- 3.5 

66.7 

CH4 

HCN 

3.8677 



3.7457 

• 3.2 

69.1 

C2H 

HCN 

3.1729 

1.1874 


3.1381 

- 1.1 

106.1 

HCN 



1.2009 





AVERAGE ABSOLUTE ERROR 




3.7 

50.9 
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TABLE IV.- CORRELATION OF BINARY DIFFUSION COEFFICIENTS FOR AN 
OXYGEN-NITROGEN-CARBON SYSTEM BASED ON DATA OF SCALA AND GILBERT 
(Ref. 3). TEMPERATURE - 12,000°R, PRESSURE = 1 ATM. 

SPECIES ^ij FROM F. -^j FROM PRES- ERROR USING ERROR IF ALL 

KINETIC THEORY FNT CORRELATION PRESENT ARE AS- 

2 2 CORRELATION SIJRED EQUAL 

i j (ft /sec) x 100 (ft /sec) x 100 (PERCENT) (PERCENT) 


f) 

02 

3. 1562 

1.0707 

3.4987 

10.8 

8.9 

Q 

N 

A. 6732 


3.4962 

-25.2 

- 26.9 

0 

M2 

3.0323 


3.4086 

12.4 

12.7 

D 

CO 

3. 1813 


3.5154 

10.5 

7.4 

0 

C02 

2.7 04 1 


2.8691 

6. 1 

26.5 

0 

C 

A. 0859 


3. 1361 

-23.2 

- 16.9 

G 

C3 

3.0408 


3.2811 

7.9 

12.5 

0 

CN 

3.3160 


3.6631 

10.5 

3.0 

02 

N 

3.3859 

1.0000 

3.7433 

10.6 

0.9 

02 

N2 

3.9602 


3.6495 

- 7.8 

- 13.8 

02 

CO 

4.0369 


3.7639 

- 6.8 

- 15.5 

02 

C02 

3. 1653 


3.0719 

- 2.9 

7.9 

02 

C 

3.0433 


3.3578 

10.3 

12.3 

02 

C3 

3.6900 


3.5130 

- 4.8 

- 7.5 

02 

CM 

4.2103 


3.9221 

- 6.8 

- 19.0 

N 

N2 

3.2441 

1.0007 

3.6469 

12.4 

5.3 

N 

CO 

3.4055 


3.7612 

10.4 

0.3 

N 

C02 

2.9097 


3.0697 

5.5 

17.4 

N 

C 

4.3167 


3.3554 

-22.3 

- 17.6 

N 

C3 

3.2668 


3.5105 

7.5 

4.6 

N 

CN 

3.5468 


3.9193 

10.5 

- 3.7 

N2 

CO 

3.9652 

1.0264 

3.6670 

- 7.5 

- 13.9 

N2 

C02 

3.1396 


2.9928 

- 4.7 

8.8 

N 2 

C 

2.9238 


3.2713 

11.9 

16.8 

N2 

C 3 

3.6396 


3.4226 

- 6.0 

- 6.2 

N2 

CN 

4.127-4 


3.8211 

- 7.4 

- 17.3 

CO 

C02 

3.1932 

0.9953 

3.0866 

- 3.3 

7.1 

CO 

C 

3.0579 


3.3738 

10.3 

11.8 

CO 

C 3 

3.7092 


3.5298 

- 4.8 

- 7.9 

CO 

CN 

4.2087 


3.9408 

- 6.4 

- 18.9 

C02 

C 

2.6492 

1.2194 

2.7536 

3.9 

29.0 

C02 

C3 

2.8901 


2.8809 

- 3.2 

18.2 

C02 

CN 

3.3357 


3.2163 

- 3.6 

2.4 

C 

C 3 

2.9458 

1.1156 

3. 1490 

6.9 

15.9 

C 

CN 

3.1740 


3.5156 

10.8 

7.6 

C 3 

CN 

3.8720 

1.0663 

3.6782 

- 5.0 

- 11.8 

CN 



0.9551 




AVERAGE 

ABSOLUTE ERROR 



8.9 

12.1 
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rigid sphere model was used when one or two monatomic species were involved. 

This procedure, in effect, means the use of different values of e^/k for 
polyatomic species, depending upon whether they are interacting with other 
polyatomic species or with monatomic species (e./k = 2670 in the former case 
in order to yield = 1-0 for the temperature considered, and e^/k 100 

in the latter case) . In order to assess the validity of this model. Table V 
lists j&. . for the hydrogen-oxygen system utilizing the values suggested by 

9^0 

Svehla, utilizing the procedure followed in Ref. 3, and using the Lennard- 
Jones model. The ^ computed using the method of Ref. 3 can be seen to 
differ by a factor of approximately two from the compilation of Svehla for 
each case involving a monatomic species . Thus one can conclude that the re- 
sults presented in Table IV represent an unrealistic test of the approximation 
procedure and should not be considered in the evaluation of the accuracy of 
the technique. It also suggests that the correlation equation (Eq. (5)) may 
have some merit for assessing the authenticity of diffusion coefficient data. 


In the preceding paragraphs correlations have been presented for four 
chemical systems. The representation of data for these systems can prob- 

ably be considered satisfactory for most purposes, in view of the uncertainties 
associated with the estimation of diffusion coefficients. The question natu- 
rally arises as to the general applicability of the correlation Equation (5). 
Further insight into this matter can be obtained by rearranging the kinetic- 
theory expression for (Eq. (4)) as follows 


"T 3 / S n.n . VK. + 

.. 1 J 1 J 

1 


1 

p a . , 3 Q. . 

LF i] i] J 

n .ft.** 
L 1 1 J 


wJ 


(17) 


where n^ can be selected arbitrarily in order to achieve a good correlation. 
As mentioned previously, the integral is unity for the rigid-sphere 

approximation. It is shown in Fig. 1 for the Lennard- Jones model, together 
with a logarithmic curve fit which is seen to be quite accurate over the range 
of T* . of interest in high-temperature boundary-layer studies 


« 1.07 (T* . ) * 1 5 9 

id id 


4 < T ij < 400 


(18) 


The 1 ^ is thus given by 


■hj 


- T 1 .6 + .159* run., (^ + 


[■(eiA)-- 0795 *" 


"(e .A)'* 0796 *" 

Pftu + Oj)* J 


"A 15 J 


L ntf J 


(19) 
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TABLE V.- BINARY DIFFUSION COEFFICIENTS FOR A HYDROGEN- 
OXYGEN SYSTEM AS PREDICTED BY THREE DIFFERENT METHODS. 
TEMPERATURE = 12,000° R, PRESSURE - 1 ATM. 


Species 

Diffusion coefficient, 
(ft 2 /sec) 

i 

j 

Recent 
Compilation 
of Svehla 
(Ref. 9) 

Lennard-Jones 
Potential Using 
Force Constants 
of Svehla 
(Ref. 8.) 

Method of 
Scala and 
Gilbert 
(Ref. 3) 

H 

H 2 

67.6 

36.0 

17.5 

H 

h 2 ° 

28.3 

26.0 

15.6 

H 

0 

27.7 

26.6 

14.1 

H 

°2 

24.6 

22.8 

11.3 

H 

OH 

29.6 

26.4 

13.2 

H 2 

H 2° 

19.6 

17.4 

17.4 

h 2 

0 

23.6 

17.7 

9.6 

H 2 

°2 

17.2 

15.0 

15.0 

CN 

X 

OH 

20.2 

17.6 

17.6 

H 2° 

0 

8.3 

7.1 

4.7 

h 2 ° 

°2 

5.7 

5.3 

5.2 

H 2° 

OH 

7.1 

6.9 

6.9 

0 

°2 

6.8 

5.6 

3.1 

0 

OH 

8.6 

7.3 

4.0 

CM 

o 

OH 

5 . 6 

5.5 

5.3 


15 
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Figure 1.- Collision integral for Lennard-Jones potential 





where k is zero for the rigid-sphere model and unity for the Lennard-Jones 
model. Equation (19) has been arranged so that the terms in brackets can 
be considered like D, 1/F^, and 1/F ^ , respectively. Note that all tem- 
perature and pressure dependence is contained in the D if the n. are con- 
strained to be constants. A successful correlation thus requires finding a 
set of n^ such that 


G(i, j) 


~n i n i (7n ± + 

(a . + a .) 3 

1 y 


h n 


( 20 ) 


is nearly constant for all species i and j in the chemical system of interest. 
It can be seen that success of this correlation cannot be predicted in general, 
since the constancy of G(i,j) depends upon the specific chemical system of 
interest. The closest thing to a general validity check, then, is to study a 
system of species selected by a random process encompassing widely different 
types and sizes of species. This has been done for a 9-component system 
(36 j) utilizing the rigid-sphere model.* The species considered and the 
resulting values of n^, chosen by trial and error, are shown in the follow- 
ing tabulation. 


i 

*L 

a . 

i 

First 

Iteration 

Second 

Iteration 

Third 

Iteration 

H 2 

2.016 

2.915 

1.0 

1.01 

1.006 

Kr 

83.8 

3.61 

1.0 

0.86 

0.820 

CO 

28.01 

3.706 

1.0 

1.01 

1.062 

CC1 4 

153.858 

5.881 

1.0 

1.08 

1.107 

SJ? 6 

146.066 

5.51 

1.0 

1.04 

1.040 

Br 2 

159.832 

4.268 

1.0 

0.80 

0.766 

AsH 3 

77.934 

4.06 

1.0 

0.94 

0.928 

C 2 H 4 

28.032 

4.232 

1.0 

1.19 

1.212 

C 2 H 5 OH 

46.048 

4.455 

1.0 

1.15 

1.165 


*The species were selected from a table of a ± from Ref. 2, pp. 1110-1112. 


17 


First, the ratio of + 77 (^ ) 2 /(o ^ + ct^) 2 was computed for each of the 36 

combinations of i and j (i.e., n^ = = 1.0). Values of were 

then assigned to each species in an attempt to improve the correlation. The 
iterative process was then repeated a second time. At this point the 36 
values of n^n^ + a j) 2 varied from 0.136 2 to 0.1816 with an 

average absolute error of 6.2 percent. This correlation is considered satis- 
factory considering the severity of the test. Furthermore, the accuracies 
could probably be improved somewhat by use of a systematic least-squares pro- 
cedure . 

Equation (17) demonstrates that the F^ can be considered as temperature 
independent if the n. are constrained to be a set of constants for a given 
chemical system. One might be able to improve correlation by allowing 
the n. and thus the F. to be temperature dependent. In order to investi- 
gate this possibility, calculations were made at several temperatures for 
diffusion coefficients based on the Lennard-Jones model. The results are pre- 
sented in Table VI in the form of an average F^ for each species and the 
percent deviation from these values. The largest single error listed in this 
table in assuming F^ independent of temperature is 0.7 percent. The average 
absolute error is only 0.13 percent, approximately one order of magnitude 
lower than the maximum errors introduced by the original approximation. Thus 
it is consistent, at least when the Lennard-Jones model is employed, to con- 
sider the F. as temperature independent and to absorb this additional error 
into the original approximation. 

As one final consideration, it can be seen from the data that have been 
presented and from Eq. (19) that the F^ tend to increase as the 77\^ is in- 
creased. A correlation of the F. obtained for the various chemical systems 
considered in this report (0-H, 0-N-C and O-N-C-H) suggested the following cor- 
relation between F^ and 77[^i 

P ± - (\)°- 461 ( 21 ) 

This correlation is shown in Fig. 2. Such an approximation is not recommended 
for general use since the F^ can be accurately determined quite simply by a 
least-squares correlation of diffusion coefficient data. The correlation 
Eq. (21) together with Eq. (5) may be of interest, however, for crude esti- 
mation of diffusion coefficients when they are unavailable. 

In summary, an approximation for binary diffusion coefficients has been 
introduced where (v - 1) v/2 . are replaced by v diffusion factors . 
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TABLE VI.- VARIATION WITH TEMPERATURE OF DIFFUSION FACTORS 
FOR AN OXYGEN-NITROGEN-CARBON SYSTEM BASED ON LENNARD-JONES 
POTENTIAL USING FORCE DATA OF SVEHLA (Ref. 8) . 


Temperature , °R 

4000 

8000 

12,000 

16,000 

Pressure , atm 

1 

1 

1 

1 

D X10 2 , ft 2 /sec 

0.6657 

2.0980 

4.1160 

6.6462 

Percent error in 
F i for: 





Fq = 0 . 7399 

-0.06 

in 

o 

0 

1 

0.04 

0.05 

F n = 1.0000 
U 2 

- .07 

- .05 

.05 

.06 

F n = 0.7378 

.04 

.03 

.11 

- .19 

f n 2 “ 1 - 0315 

.04 

.03 

.12 

- .19 

F co = 1.0217 

- .04 

0 

.07 

- .03 

f co 2 “ i- 2690 

- .19 

- .18 

.12 

.25 

F c = 0.6692 

.42 

.35 

- .70 

- .06 

F^ = 1.0928 
c 3 

- .15 

- .16 

.04 

.27 

f cn " 1 - 0343 

.03 

o 

* 

.12 

- .17 
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Diffusion Factor 



Molecular Weight,#^ 


Figure 2. - Correlation of Diffusion Factors with Molecular Weight. 
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The primary advantage of this approximation is that the diffusional mass flux 
of species i can then be expressed explicitly in terms of gradients and prop- 
erties of species i and of the system as a whole. This, in turn, permits 
substantial simplification of the boundary-layer equations. The F i are 
independent of pressure and are only slightly temperature dependent. Further- 
more, they can be considered to be independent of concentrations. This has 
the advantage that the F^ can be determined a priori for a given set of 
molecular species, becoming simple input data in a numerical solution. It 
was shown that the actual binary diffusion coefficients are approximated 
within a few percent with the following exception. When the molecular weights 
of species i and j differ by an order of magnitude, the 3 . . determined 
by the least-squares curve fit can differ from the actual . by as much as 
15 percent or so. It can thus be concluded that the approximation procedure 
is quite accurate for the majority of diffusing pairs and is within the range 
of certainty of diffusion coefficients for most if not all other species. 


SECTION 3 
MIXTURE VISCOSITY 


The viscosity of a multicomponent mixture according to rigorous first 

2 

order kinetic theory is given by Hirschfelder , Curtiss, and Bird as 


H 11 H 12 ' 


H 12 H 22 


H lv H 2v ' 


H lv H 2V 


. H 


lv 


. . . H 


2v 


vv 


X 2 

* * 

* X v ' 

'll 

H 12 ’ 

• • H lv 

'l2 

H 22 * 

• * H 2 v 

* 

* 

• 


. H 


vv 


( 22 ) 
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where x i is the mole fraction of species i, v is the number of molecular 
species, and the and are given by 


H. . 

IX 


= ^ + 7 


2x i x k 


(% + 




RT 

P ^ik 


k/i 


3 ^k 

+ — — A* 

5 7n ± ik 


(23) 


H ij 


2x . x . 

1 J , .-. 

W; + in 


RT 


~ r 1 - -i . i 

ij L 5 i 3 J 


(i ? j) 


(24) 


with P the pressure, T the temperature, R the universal gas constant, 
and the molecular weight of species i. is the binary diffusion 

coefficient for pure species i and j given by Eq. (4) , is the viscosity 
of the pure species i which can be expressed in terms of the self-diffusion 
coefficient as 


A *xj 


is defined as 


6A ii P ^ i;L 


A *ij 




IJ 


(25) 


(26) 


p is the density, T* ^ is a reduced temperature given by Eq, (15) , and 

of 1 . #1 ^ * (T* .) and oj* ' S ^ * (T^j) are integrals for calculating the transport 
coefficients (given for several molecular potential models in Appendices to 
Ref. 2, for example) . 

Hirschfelder, et al 2 show that Eq. (22) can be expanded as 


^mix 



i=l 


U V 

EE 

i=i j=i 


X i X 1 H ij 


+ 


u u u 



i=l j=l k=l 
j^i k^i 


x i x k H ii H ix 
H ii H 


:A* 


(27) 
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The primary contribution to the mixture viscosity is given by the first term 
in this expansion which represents the diagonal terms in the determinants of 
Eq. (22) . In fact, all but the first term vanish for = 5/3. With this 
approximation, the expression for mixture viscosity becomes 


U • 
^mix 


I 




Z 

j=i 


2x . x . 


RT 


i j 


(28) 


This equation does not predict mixture viscosities well because of the unrea- 
listic value* assumed for h*. . and, consequently, the neglect of higher 

x 3 n 

terms in the expansion of Eq. (27) . Buddenberg and Wilke have shown by an 
extensive analysis of experimental data that Eq. (28) describes mixture vis- 
cosity if the theoretically determined constant "2" is replaced by an empiri- 

2 

cal value "1.385". Consequently, Hirschfelder et al recommend that a good 
approximation to Eq. (22) is given by the semi-empirical equation 


v r 


mix 


-I 


i=l 


X i^i 


R ^i ~ 

x, + 1.385 -—± V 
i P ty. *-> 


V 


1 j=l 

j^i 


(29) 


Introducing the correlation equation for diffusion coefficients (Eq. (5)) 
into Eq. (29) yields 


^mix 


V 

I 


i=l 


X i^i 


x. + 1 
i 


.385 — 


Vi 


— m. 

PD l i 


j=l 

j^i 


x . F . 
3 3 


(30) 


Utilizing the definitions of the system properties (i^ and \jl ^ given by 
Eqs . (9) and (10), respectively, Eq. (30) becomes 


* A better value for 


A* . 
13 


is 1.1 to 1.2. 
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( 31 ) 


mix 


V 

-I 


X .|1 . 
11 


i=l 


RT 

x . + 1.385 ( Ml _ x .f.) 


The equation of state for the mixture is given by 


P 


R 

% 


p T 


(32) 


where W[ is the molecular weight of the mixture. Substituting Eqs. (32) f 
(25) and (5) into Eq. (31) and rearranging yields the following result for 
the viscosity of the multicomponent mixture 


mix 


_ s»y 
^i l . 


x i*i 


i=l 


1.385 + 


x i F i 


6A* . 

ii 




1.385 


(33) 


It should be noted that calculation of the mixture viscosity thus has been 
reduced to a single summation involving system properties which appear 
elsewhere in the boundary-layer equations 


SECTION 4 

MIXTURE THERMAL CONDUCTIVITY 


The thermal conductivity in a polyatomic gas mixture has been shown by 
12 

Hirschfelder to be represented by 

\ . = \ . + . (34) 

mix mono-mix int 

where \ . is the thermal conductivity in a mixture computed neglect- 

mono-mix J 

ing all internal degrees of freedom and is the contribution to the 

thermal conductivity of the mixture due to the internal degrees of freedom 
of the molecules. 
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The \ . is given by Hirschf ".der et al as 

mono-mix 


X = 4 

mono-mix 


'll 

L 12 

• * * L lv 

X 1 

J 21 

L 22 

• • • ^2v 

x 2 

• 

J vl 

L v2 

• • * L vv 

♦ 

• 

X 

V 

: 1 

X 2 

. . . x v 

0 


L 11 L 12 
L 21 L 22 


• • L i 


• • L , 


where 


L ii = 


L vl L v 2 • • • L w 


16 T y x i x k [•¥• V + jr V - 3 V B !k + 4 WlJ 

25 p ^ + v^ik 


„ 16 T 

L ij - TTT 


x.x. &!.&!• 55 -| 

£ 1 3 3 ¥■ - - 4A* (i / j) (37) 

P (#!; + L 4 13 13 J 


i J ID 


with X imono the tlier][nal conductivity of the pure species neglecting all in- 
ternal degrees of freedom of the molecules 


n j_ u 

4 


and defined as 


5 n iV 8 ) *( T h) - 4 fi iV a) * (T !i ) 

iJPV) 

ii v ij 
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Approximate expressions can be developed for thermal conductivity in a 

manner analogous to the procedure for viscosity described in Section 3. 

2 

Several formulations are referenced by Hirschfelder. For the present pur- 

13 14 

pose a method suggested by Mason and Saxena and discussed by Brokaw is 
adopted. However, the same procedures could be applied to other rela- 
tions. 

Mason and Saxena showed that an equation analogous to Eq. (29) was 

applicable for X . by replacing the |j. which appears in the numer- 

mono-mix i 

ator by ^i m0 no and replacing the empirical coefficient 1.385 by 1.385 x 
1.065, that is. 


mono-mix 


-I 


i=l 


Vi 


mono 


x^ + 1.065 x 1.385 


RTp . v x. 


j-1 

j/i 


(40) 


Substituting for Xi „ 
x mono 

from Eq. (5) yields 


from Eq. (38) , for from Eq. (25) and for 


mono-mix 


V 

-E 


i=l 


6A ? . 4 


15 - R x i 

— pD W L Fp 


x. + 1.065 x 1.385 + -r|=- rr=- Y x.F. 

l p 6^ F i ^! i Lj 3 3 


j-1 

j^i 


(41) 


Using the equation of state (Eq. (32)) and the definition of \i ^ (Eq. (9)) 
and rearranging yields 


mono-mix 


v 

15 RfiD V" 
4 


i=l 


X i /F i 


1.065 x 1.385 + 


x i F i 


I 6Af, % 


n " 1-065 


x 1.385 


(42) 

The contribution to the thermal conductivity from the internal degrees 

12 

of freedom X . is given by Hirschfelder as 
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X. ^ 

int 


V 

-I 


mono 


i=l 


v 

£ 

j=i 




(43) 


where ^i mono i s given by Eq. (38) and X ^ is the true thermal conductivity 

of the i^rx species. An expression for X ^ taking into account the internal 

structure of a polyatomic molecule was first given by Eucken. An improved 

2 

relation suggested by Hirschfelder et al has the form 


Xf - 


1 

4 


15-6 


P Ai 


- 15-10 


p .3. . 

H i n 


v i 


(44) 


where is the density, Cp^ the specific heat at constant pressure, and 
C V £ the specific heat at constant volume for the pure species i. Utilizing 
the relations 


and 




R 

% 



it follows that* 




h. 


mono 


\ 

p ^ii 7tl 


1 R 
2 


% 


(45) 


(46) 


(47) 


Substituting into Eq. (43) and utilizing the correlation Eq. (5) yields 



(48) 


♦This form is suggested directly by Hirschfelder as Eq. (59) of Ref. 12. 
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It is shown in Ref. 6 that in addition to and u ^ f it is convenient 

in multicomponent boundary- layer problems to define 


Y Z.C 

zd 1 Pi 


( 49 ) 


i-1 


and 


V" fi _ ir!i 

L \ ^ L f i 

i=l i=l 


(50) 


where is defined by Eq. (3) . Substituting these into Eq. (48) yields 


X . 


pDd, 


"int ) | C p ~ 2 ^3] 

The mixture conductivity is thus given by the sum of Eqs. (42) and (51) 


(51) 


X ^ PP 

mix | < 


f - 

r X i 15 R 



Jy 

f. 4 m 



\k 

x . F . . 

1.065 x 1.385 + | 

M 1 

6A* . 
li 

5 

- 1.065 x 1 . 385 | 


^2 ~ 5 

+ W C p “ 2 ^3 


(52) 


As with mixture viscosity, the mixture conductivity can thus be computed by 
a single summation involving no new system properties other than which appear 
elsewhere in the bound ary- layer equations. 
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SECTION 5 

MULTICOMPONENT THERMAL DIFFUSION COEFFICIENT 


The multicomponent thermal diffusion coefficient is given by Hirschfelder , 
2 

Curtiss, and Bird as 


D. 

i 


8 ^. 

5R“ 
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T oo 
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T 01 
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00 
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(53) 


where 


L 00 


= o 


(54) 
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3 ¥■ k 


(55) 



8 T 
5 P 


v pva ll c i« - i ) ~ 

L L <V°V 3 j7 _ 

X =1 


(56) 



8 T 
5 p 


gi iijik - 1 
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(57) 


L 


10 

jk 


\ t 01 
VI, L jk 


(58) 


L 11 and lP; (j ^ k) are the same as L. . and 
3D D K 

the equation for thermal conductivity of a mixture of 
(36) and (37)) with a change of subscript from i to 
k to X , and 


which appear in 
monatomic gases (Eqs. 
j, j to k, and 


c h 


fi* 1 ' 1 )* 


< T ti> 


(59) 


This expression could be simplified somewhat by employing the approximation 
for (Eq. (5)) into Eqs. (36), (37), and (55) through (57). To illus- 
trate, Eq. (55) can be reduced to 


.00 

J jk 


16 T 
25 PD 


F -i 

< j F j x k F k + x k\ %- 2^ 


x x F i 


i=i 


Introducing |i^, defined by Eq. (9) , yields 


00 _ 16 T 

* = 25 Pi 


x j F j x k F k " \ ^ ) + M l x k\ ^ 


(60) 


(61) 
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However, the calculation of multicomponent thermal diffusion coefficients 

would still be algebraically quite complex. Therefore, it would be very 

T 

desirable to develop an approximate relation for the , if possible. One 
might expect at least fair success since the situation is somewhat analogous 
to the problem of ordinary (concentration) diffusion. As in the case of the 
2Kj, the is independent of the diffusive and thermal fluxes. in the 

former case it was found (see Eq. (5)) that one binary diffusion coefficient 
can be used to represent unequal diffusion coefficients for all species with 
the unequal diffusion effects taken into account through the diffusion fac- 
tors F.. Thus it might be possible to represent multicomponent thermal dif- 

1 — t 

fusion by a single binary thermal diffusion coefficient D with the effects 

of the unequal thermal diffusion coefficients taken into account through 

thermal diffusion factors analogous to the . If one is fortunate, 

in fact, it might turn out that the are the same as the F^ or closely 

related, but this is not a priori evident nor is it necessary for a successful 

approximate formulation. 

On a microscopic basis, the mean thermal diffusion velocity of a statis- 

T 

tical set of molecules, , must result from some characteristic of this 

set as it relates to the remainder of the system. Since this set can be con- 
sidered extremely small on a macroscopic basis, the "remainder of the system" 
can be considered as the system as a whole. Thus one can postulate 

V ± T * (G i ■ G > (62) 

where G represents some mean value of the property G^ for the system as 
a whole. The thermal diffusion coefficient represents a total mass flux of 
species i, therefore 


°i T “ p K i (G i ~ G) (63) 

where pK^ is the mass concentration of species i. It is required, however, 
that 


u 

£ D ± T = 0 (64) 

i=l 

Equations (63) and (64) lead to the requirement 

G - Y. K j S j 1651 

j=l 
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Thus the following relation is proposed as a basis for correlating thermal 
diffusion data 



( 66 ) 


where, as mentioned earlier, D is a property of the system as a whole and 
is a thermal diffusion factor for molecular species i, analogous, respec- 
tively, to the D and the for ordinary diffusion (see Eq. (5)). 

In order to evaluate the validity of this expression and, in particular, 
to ascertain appropriate expressions for and D , recourse has been 

taken to thermal diffusion data and calculations available in the literature. 
As a logical first step in the verification of the proposed model, data for 
binary systems have been investigated. 

Thermal diffusion data for binary systems are sometimes presented (e.g.. 
Ref. 2) in terms of a thermal diffusion ratio k T defined as 


T 

k a 

T W P-»ab 


(67) 


The quantity k T is so defined such that when k T > 0 component a moves 
to the cold region, and when k T < 0 component a moves to the hot region. 
Consistent with the usual convention species a is herein considered to 
be the heavy species so that under normal conditions k^ is usually positive. 

A second and more useful form for the presentation of binary thermal dif- 
fusion data is the thermal diffusion factor, a^' defined by 


a 


ab 



K aV^ab 


( 68 ) 


This has the advantage over k that it is much less concentration dependent. 

r T 

In terms of a , , the binary thermal diffusion coefficient D is expressed 
ab u 

as 


D a T = a abP K aVab 


(69) 
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The can be eliminated from this equation in favor of D and 
the use of Eq. (5) , which is exact for a binary gas mixture if the 
properly chosen. This yields the result 


^ K a K b 


“ a abP D F7FT 

a b 


F i 


by 

are 


(70) 


Likewise, the correlation equation proposed in Eq. (66) becomes for a binary 
system 


°a T = dT ’ K aV G a " G b> ( 71 > 

or from Eq. (70) 

a ab = =T (G a - ^ F a F b < 72 > 


Use will be made of these relations later. First, however, it is pertinent 
to review the behavior of ot^. 

There are several available theoretical approximations for a , . Mason 15 

16 

suggests that the method of Kihara be used in preference to the other avail- 
able approximations. The Kihara equation can be written in the form 


where 
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with 


a 




D* = 




-Q 


ab 


< T 3b> 


a \ s 

aa 1 


ab , 


(77) 


and C i j ' and a ij tlie same as defined previously. The relations 

for Q^, and D£ are obtained from the above by interchange of sub- 

scripts. As pointed out by Mason, the principal contribution to the tem- 
perature dependence of comes from the factor (60*^ - 5) , whereas the 

term involving the S and Q determines the composition dependence of a T , 
but depends only slightly on temperature. This is an important feature which 
will be utilized later. 

Mason'*’ 5 presents computed using Eq. (73) for several binary mix- 

tures of gases (primarily noble gases) as; 1) a function of concentration 
(from 0 to 100 percent) for a constant temperature (327°K) (his Fig. 5) and 
2) a function of temperature (0 to 600°K) for equivalent mole fractions (nomi- 
nally 50 percent) (his Figs. 3 and 4) . Typically, the following observations 
can be made: 


1. oc^k increases from a finite value at 0 percent light gas to some 
higher value at 100 percent light gas in a manner which appears to 
behave consistently with molecular weight ratio ? 

2. At a given mole fraction, a &b increases as the ratio is 

increased (recall that a refers to the heavy gas); 

3. The a ^ tends to become quite insensitive to temperature for tem- 

peratures above 300 to 800° K depending upon the particular pair of 
gases under consideration. This is consistent with the observation 
that the primary temperature dependence of is contained in 

the term (60*^ - 5) , since C*^ tends to become uniform at ele- 

vated temperatures.* 

On the basis of these observations, it is seen that a^, as a first 
approximation, can be considered to be independent of temperature for tem- 
peratures above 1000° K or so, depending only on concentration and molecular 
weight. In anticipation of the desirability of forming a relationship be- 
tween the G. and the F., correlation was attempted with the ratio F /F, 
li a d 


*For example, the C* for the Lennard-Jones potential is shown by Hirsch- 
felder et al 2 to be at> 0.9483 for T* of 30 to 400 (the highest value of 
T* tabulated) . 
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rather than since, as shown in Section 2, the diffusion factor varies 

approximately with molecular weight as 




%. \ ° ,46 
a 

Kl 


(78) 


An investigation of the boundary- layer equations revealed that tremendous 
simplification results if 



(79) 


with c^_ a constant. The binary thermal diffusion calculations of Mason 
are correlated with Eq. (79) in Fig. 3. 

The data points correspond to molecular weight variations of 1.98 to 

32.8 for the following mole fractions of the heavy species; x^ — ► 0, x ^ 0.50, 

a a 

and x^ — ► 1.0. These were obtained in the following manner. First the asymp- 
totic values of at high temperature for & 0.50 were estimated di- 
rectly from Fig. 5 of Ref. 15. Then values of C a a £> ( x _ 327 °k and 

I- a ab^ x 0)] T = 327°K were obtained from Figs. 3 and 4 of the same refer- 

ence. These were then corrected in an approximate manner by multiplying by 
the ratio of |> ab (x *. 0.50) ] Asympt> to [a gb (x « 0.50) ] T = 327 o K . These 
calculations are summarized in Table VII. 

It can be seen that the data for the absolute value of the abscissa less 
than unity correlates quite well with a straight line through the origin* with 
a slope of -0.5 corresponding to c^ = -0.5. The scatter is of the same order 
that experimental data shown by Mason scatters from his calculated values. 

The correlation appears to begin to fail for values of (1 - F /F, ) /I” 1 - x (1 - 
WH of unity, overpredicting a ^ by about 30 percent for the three 
data points shown in this range. In order to evaluate the seriousness of 
the relatively poor correlation at large values of this parameter, the param- 
eter is presented in Fig. 4 in terms of and x a - Recalling that the 

subscript a refers to the heavy species, it can be seen that the thermal 
diffusion of a small concentration of a species with large molecular mass 


*It must go through the origin as a , must vanish when molecular weights 
are equal. a 
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Data of Mason (Ref. 15) 
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Figure 3. - Correlation of Binary Thermal Diffusion Data for G.= 1/F . 





TABLE VII.- CORRELATION OF BINARY THERMAL 
DIFFUSION DATE OF MASON (Ref. 15 ) . 
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Figure 4.- Variation of binary thermal diffusion correlation parameter 
with concentration and molecular weight ratio 



relative to the gas mixture (say x & < 20 percent and > 1°) or a 

large concentration of a species with low molecular mass relative to the 
gas mixture are not accurately represented. Fortunately, this is not a 
serious shortcoming for most boundary-layer problems. The former is not a 
practical problem because few boundary- layer species have molecular weights 
10 times greater than the system molecular weight. The latter is unrealis- 
tic, since a large mole fraction of a low molecular weight species would 
not permit a large molecular weight ratio betwen that species and the system 
as a whole 

A study of the theoretical work of Laranjeira^ 7 (see Appendix C) leads 

to the conclusion that the a , should behave as 

ab 



The of Mason are correlated in accordance with this equation in Fig. 5. 

The data appears to correlate better for large values of but the scatter 

at lower values is somewhat greater. On the basis of this cursory study, 

Eq. (80) appears to offer little if any improvement in the accuracy of the 
correlation. Therefore, since Eq. (79) leads to substantially greater alge- 
braic simplicity, it is tentatively recommended as the correlation equation. 
Further studies should be performed with additional thermal diffusion data. 

Comparison of the correlation Eq. (79) with the proposed relation for 
binary systems (Eq. (72)) indicates 


G a = !/P a . G b = 1/F b 


(81) 


5 T 


C t D 


x F + x, F, 
a a b b 


(82) 


A rigorous extension of these results to multicomponent systems would require 
the investigation of multicomponent thermal diffusion data. Such data could 
be generated by parametric solution of Eq. (53) . For the present, however, 
it is noted that the particular form of Eqs . (81) and (82) suggests the fol- 

lowing generalizations for multicomponent systems 
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ation of Binary Thermal Diffusion Data for G 1/F. 


and 


C t D 


y x .f . 

{-> 3 3 

j=l 


(84) 


The theoretical work of Laranjeira 17 appears to substantiate this generaliza- 
tion (see Appendix C) . 

Substituting Eqs . (83) and (84) into the proposed correlation Equation 

(66) yields the result 


D. 

l 


_ r_L _ £ Hii 

pDK. If. ■ f-i, F . 
t K i L i i=l t-I 


(85) 


E 

j=l 


. F . 
' x 


Utilizing the previously defined quantity ix^ (Eq. 9)) yields 


D. 

i 


c t pDK i 


f - Y. 

■ 1 j=i 


Hi 

F j 


(86) 


Furthermore, using the definitions of lij (Eq. (10)) and (Eq. (8)) 

yields the ridiculously simple result 


T c t pEU 2 

D i = TTIT" (z i - V 


(87) 


Upon substitution of Eq. (87) into the expression for diffusive flux of 
til 

the x molecular species, j ^ , and into the energy equation, the processes 
of thermal diffusion (Soret effect) and diffusion thermo (Dufour effect) are 
both taken into account with surprising ease, thus allowing the inclusion of 
these effects with effectively no increase in computational time compared to 
their complete neglect. To illustrate, the expression for mass diffusive 
flux for assumed independent of temperature (Eq. (12)) simplifies to 


P°^2 

I~ az i L „ ) i 

[ a in n 2 , 3 in T 1 


W + ( Z i ~ K i) 1 

I ay t a y ) 

— 


( 88 ) 


This relation is developed in Ref. 6. Modification of the boundary-layer 
equations to include this approximation for multicomponent thermal diffusion 
coefficients is also presented in Ref. 6 . 
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APPENDIX A 


DEVELOPMENT OF APPROXIMATE RELATION FOR 
MASS DIFFUSION FLUX 


Substituting the approximation for binary diffusion coefficients (Eq. 
(5)) into the Stefan-Maxwell relation (Eq. (2)) and rewriting in terms of 
mass fractions yields 



where, for convenience, a total diffusional mass flux has been defined as 
the sum of the molecular and thermal diffusional fluxes. 


= j . + D. 
J x x 


T l Un T 

dy 


(A- 2) 


Multiplying each side of Eq. (A-l) by F^, summing over all i, and 
noting that the sum of the diffusive fluxes is zero and the sum of the mass 
fractions is unity yields: 


Z 


J.F. 



— <r — 1 7H. 3x. 

= £D \ _i i 

v? L F ± 

i 


£D 


I 


Fj 9y 


Substituting Eq. (A-3) into Eq. (A-l) results in 


(A- 3) 


dx 


F i y 

\ L, F j ay 


i . SL- i?A V -ifi 

“ = PD *i L 


(A-4) 


At this point it is convenient to define several new quantities. 


Z. H ^ i x./F.n 2 


(A-5) 


'1 ' Z X i F J 


(A-6) 
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(A-7) 


£ ^j x j /F j 


I*. 


</Fj S ) (dPj/dT) 


(A-8) 


Multiplying Eq. (A-5) by p 2 and differentiating the result with respect 
to y yields 


\ 3x i ^ iXi 3F i 

F~ Sy F 8 Sy" 


az. 

m 2 ay 1 + z i aF 


(A-9) 


Substituting Eqs. (A- 2) and (A-5) through (A-9) into Eq. (A-4) and noting 

that V Z. =1.0 and that the F. are functions of temperature only 
*-* D 1 

j 


h + 


i 

T 


aT pd 

ay m. 


^2 

az. 

i , 

(Z. - K.) 
' 1 1 

±. 4 _ t r 

f-L. ffi . a \ 

aT 


ay 


_ T • 

ay i 

| P .» « 4 

ay 


(A-10) 


This is the desired explicit relation for j ^ in terms of gradients and 
properties of species i and of the system as a whole (Eq. (7) of report) . 
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APPENDIX B 


A FILM-COEFFICIENT APPROXIMATION FOR 
MULTICOMPONENT DIFFUSIONAL MASS TRANSFER 

In this Appendix, an approximate relation for evaluating diffusional mass 
transfer rates in a multicomponent gas mixture is presented, the magnitude of 
the errors introduced by the approximation are assessed, and a film-coefficient 
relation for extending the expression to mass transfer in a multicomponent chem- 
ically reacting boundary layer is proposed. The film-coefficient approximation 
is based upon arguments introduced in Ref. B-l, and represents an extension to 
the results presented therein. The approximate relation for mass diffusion 
flux in a multicomponent mixture is presented in Section B.l. A description 
of the proposed film-coefficient approach is presented in Section B.2. 


B.l A MODIFIED FICK'S LAW FOR DIFFUSIONAL MASS FLUX IN A MULTICOMPONENT 
MIXTURE 


It was shown in Appendix A that the bifurcation approximation for binary 
diffusion coefficients yields the following explicit relation for the diffu- 
sional mass flux in a multicomponent gas mixture, in the absence of thermal 

and pressure diffusion and if the diffusion factors, F^, are assumed indepen- 
dent of temperature; 


pDu 2 pZ. (Z i - K i ) cHi 2 " 
dy + ~ 2 dy - 


(B-l) 


In the event that all binary diffusion coefficients, ^ , are equal (J^ = ^ 2 ) , 
the F ± can be set equal to unity. Then Z ± = K ± , D = J& 12 , U 2 /d x - W\, and 
Eq. (B-l) reduces to Fick’s Law: 


a Ki 

= - ^12 sr (B_2) 

It is noted that the form of Fick's law provides the motivation for express- 
ing the diffusional flux in terms of a film coefficient (involving the product 
p^l 2 ) and a driving potential, AK^ . 

Examination of Eq. (B-l) for multicomponent diffusional mass transfer 
reveals that an equation similar to Fick's law (Eq. (B-2)) results if the 
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second term in Eq. (B-l) can be ignored. The relative magnitudes of the 
first and second terms in Eq. (B-l) have been examined for several multi- 
component gas mixtures and it has been found that the second term is indeed 
relatively small. The accuracy checks were performed in the following manner. 

The first step was to select a chemical system, mole fractions, x^ , and 

diffusion fluxes, j., satisfying the requirements that S x . =1 and 
c"* 1 1 1 

4j j ■ =0. The mole fraction gradients, dx./dy, were then computed from 
i 1 i 

the Stefan-Maxwell relations (neglecting thermal diffusion) 







Recalling the definitions. 


Z. 

i 


Vi 

F i^2 


and 



j 


(B-3) 


(B-4) 


(B-5) 


enabled forming an expression for the "Z" gradient in terms of the mole frac- 
tion gradient. 

Vi 
ay F i d 2 

' j 


_1_ _i JL \ 3 

x. 5y 


M 2 Fj By 


(B-6) 


The 11 Z" gradient was then computed from Eq. (B-6), employing the specified 
mole fractions and the mole fraction gradients computed from Eq. (B-3) . The 
F^’s were taken to be proportional to ^. 2 in these studies. Neglecting the 
dfi 2 /dy term in Eq. (B-l) , the diffusional mass flux is given by 


] i 


pD U 2 3Z i 

Wiip ay 


(B-7) 
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The approximate j ^ computed by Eq. (B-7) were then compared to the spe- 
cified to assess the accuracy of Eq. (B-7) . 

Several chemical systems were examined as described above. The results 
of one of the more severe tests are presented in Table B-l for the system 
0, CO, and C0 2 (representing a molecular weight variation of 2 to 44) . 
Based upon the average absolute error of 6.85 percent and maximum error of 
18.37 percent, it can be concluded that Eq. (B-7) is reasonably good for 
representing diffusional mass transfer in multicomponent gas mixtures. It 
is noted that when all diffusion coefficients are assumed equal, Eq. (B-7) 
reduces identically to Fick's law (Eq. (B-2)). 

It is emphasized that Eq. (B-l) is more accurate than Eq. (B-7) and, as 
such, should be employed for detailed computation of mass transfer events. 
When an approximate treatment is adequate, such as for many engineering ap- 
plications, an integrated form of Eq. (B-7) may be employed for detailed 
computation of mass transfer events. When an approximate treatment is ade- 
quate, such as for many engineering applications, an integrated form of Eq. 
(B-7) may be employed with fair accuracy. The integrated form or "film 
coefficient" representation of diffusional mass transfer processes in the 
multicomponent boundary layer is considered in the next section. 


B . 2 A FILM-COEFFICIENT APPROACH 

In this section, the approximate relation (Eq. (B-7)) developed above is 
utilized to formulate a film-coefficient representation of mass transfer in 
the multicomponent boundary layer. An "Ohm's Law" type expression of the 
form, j ^ = CAP is desired, where C is a film-coefficient and AP is a 
suitable driving potential for mass transfer. A form for the driving poten- 
tial is rationalized in Section B.2.1 whereas an approach for evaluating 
the film coefficient is presented in Section B.2.2. 

B.2.1 The Driving Potential 

In order to rationalize a form for the driving potential it is appro- 
priate to consider the boundary- layer conservation equation for molecular 
species in the absence of thermal diffusion. 

dK. dK. / dZ. \ 

P u ai i+ P v aF = fy |p*eff ¥T J + *i (B_8) 
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TABLE B-l 



ERRORS IN 

DIFFUS IONAL 

MASS FLUX INTRODUCED BY 



NEGLECT 

OF b[i 2 /by 

TERM IN EQUATION 

B-l 


Species , 
i 

Mole Fraction, 

x . 
i 

Diffusion 
Flux, j. 

Diffusion Flux, 
Neglecting 
byi^/by Term 

percentage Absolute 
Error in j ^ With 
o\i 0 /by Term Neglected 

H 2 

.076923 

- .225000 

- .234988 


4.44 

O 

.076923 

-.1 . 025000 

-1.013632 


1.11 

CO 

.769230 

.375000 

.338282 


9.80 

C0 2 

.076923 

.875000 

.910338 


4.03 

H 2 

.062500 

- .225000 

- .233771 


3.89 

0 

.062500 

-1 . 025000 

-1.022300 


0.26 

CO 

. 625000 

.375000 

.306095 


18.37 

C°2 

.250000 

.875000 

.949976 


8.56 

h 2 

.062500 

- .225000 

- .233236 


3.66 

o 

. 625000 

. 375000 

. 331520 


11. 59 

CO 

.250000 

-1.025000 

- .997998 


2.63 

C°2 

.062500 

.875000 

.899714 


2.82 

H 2 

.052631 

- .225000 

- .232314 


3.25 

0 

. 526315 

.375000 

. 315384 


15.89 

CO 

.210527 

-1.025000 

-1.017707 


0.71 

co 2 

.210527 

.875000 

.934637 


6.81 

H 2 

.250000 

- .225000 

- .259263 


15.24 

o 

.625000 

.375000 

. 360266 


3.92 

CO 

.062500 

-1.025000 

-1.011923 


1.27 

C°2 

.062500 

.875000 

.910921 


4.10 

«2 

. 210527 

- .225000 

- .256334 


13.91 

0 

.526315 

. 375000 

. 322715 


13.94 

CO 

.052631 

-1.025000 

-1.020750 


0.41 

C°2 

.210527 

.875000 

.954369 


9.07 

H 2 

.210527 

- .225000 

- .254236 


12.97 

0 

.526315 

.375000 

. 346711 


7.54 

CO 

.210527 

-1.025000 

- .992857 


3. 14 

C°2 

.052631 

.875000 

. 900383 


2.89 

h 2 

. 181818 

- .225000 

- .251650 


11.82 

o 

.454546 

. 375000 

. 327074 


12.77 

CO 

. 181818 

-1.025000 

-1.013292 


L. 14 

C°2 

. 181818 

.875000 

.937868 


7.18 

Average 

Absolute Error ( 

% 



6.85 
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In this equation, diffusional transport is represented by the approximate 
Equation (B-7) and an effective diffusion coefficient has been defined as 


£ 


eff 


D ^2 

m ii. 


(B-9) 


The frozen, or "nonreacting" , case will be considered first, for which the 
creation term, =0. 

It will be useful for subsequent analogy arguments to consider also the 
boundary-layer species conservation equation in the absence of chemical reac- 
tions and with all diffusion coefficients equal. 

3K i 3K i 3 

pu _ + pv _ = — 

For simple Couette flow (with constant properties and no pressure gradient) , 
the terms on the left-hand side vanish and Eq. (B-10) may be integrated di- 
rectly to yield the species mass flux at the wall, jj_ , in terms of the 
channel height, 6. 

h = Pe U e C M (K i - K i > (B - 1L 

w we 






12 dy 


(B-10) 


where the simple Couette flow film coefficient for mass transfer, p u C w = 
p *^12^“ Thus, for simple Couette flow, the mass-transfer coefficient is 
related directly to the diffusion coefficient, &Y2* When the convective 
terms on the left-hand side of Eq. (B-10) cannot be ignored, the film-coef- 
ficient driving-potential Eq. (B-ll) is still appropriate for calculating 
the species mass flux. The driving potential is clearly represented by the 
mass fraction difference across the boundary layer - K^ e ) , but the 

film coefficient must now include the effect of mass transfer by both con- 
vective and diffusive mechanisms. The results of a large number of experi- 
mental investigations and data correlations are summarized and presented by 
Spalding (Ref. B.2) where it may be noted that for boundary-layer flows the 
mass- transfer coefficient depends approximately on other 

things remain unchanged. The actual exponent varies between 1/2 and 3/4 
for the wide range of conditions considered, the 2/3 exponent being a rep- 
resentative average value. Thus for boundary-layer flows, the diffusion 
events, represented by the right-hand side of Eq. (B-10), depend upon the 
2/3 power of the diffusion coefficient and significant mass transfer 
results from convective events, represented by the left-hand side of Eq. 
(B-10) . 
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Examination of Eq. (B-8) for nonreacting multicomponent Couette flow 
suggests a film-coefficient driving-potential relation for mass transfer 
similar to Eq. (B-ll) , but with a "Z" rather than a "K" driving potential. 


] . = p u c u (Z. 

K e e M 1 
w w 


2 i > 

e 


(B-12) 


where the mass-transfer coefficient is obtained by direct integration of 
Eq. (B-8) with the left-hand side set to zero. 


P 


u c„ 

e e M 


gfeff 

6 


(B-13) 


Equations (B-12) and (B-13) clearly represent a solution to the multicomponent 
nonreacting species conservation equations for Couette flow. Attention will 
now be directed to species conservation in the boundary layer. 

Examination of Eq. (B-8) and consideration of the above discussion sug- 
gests that dif fusional mass transfer in a multicomponent nonreacting (i|k = 0) 
boundary layer may be represented by a " Z" potential and convective mass 
transfer by a mass fraction potential. It would seem reasonable to expect 
that the overall mass transfer resulting from both diffusional and convective 
events can be represented by some kind of an average potential involving both 
Z ^ and . For example, the driving potential could be represented by 
A(Z^Y) where the exponent, y, would depend upon the relative significance 
of diffusional transport for the particular situation. For Couette flow, 
y =1.0, whereas for boundary- layer flows, y < 1.0. The proper value of y 
for boundary-layer flow is not apparent and may be established with confidence 
only after testing the postulated potential by correlating either experimental 
data or exact numerical solutions of the boundary- layer equations. Recalling 
the earlier discussion on boundary- layer mass transfer for equal diffusion 
coefficients, it is postulated that the relative contribution of mass transfer 
by diffusional and convective processes can be taken into account by use of a 
2/3 power dependence on the diffusion coefficient. On the basis of this 
observation, a y of 2/3 in the Z, K driving potential would take proper 
account for the relative mass transfer by diffusional and convective events. 
The film-coefficient equation for mass transfer in the multicomponent, non- 
reacting boundary layer may then be written as follows 

j. = p u C, (Z* - Z* ) (B-14) 

J i H e e M i l v ' 

w we 


i 
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where 


Z* 

1 



1 X 

z 3 K r v 


j 


( B— 15 ) 


and y = 2/3 for boundary-layer flow. The normalizing summation is included 
so that the Z* will sum to unity. 

The above film-coefficient approach for the nonreacting boundary layer 
may be extended to the reacting boundary layer in the following manner. By 
applying the Shvab-Zeldovich transformation to Eq. (B-8) (multiplication by 

B— 3 

the mass of element k in molecule i, , and summation over all species 

the creation term, i|k # is eliminated. As a consequence, the resulting equa- 
tion for conservation of chemical elements in a reacting boundary layer takes 
the same form as the molecular species conservation equation for the nonreact- 
ing (frozen) boundary layer 


where 


dK k Sit az* 

pU as - + pv IT" = 3y p ^eff ly - 




(B-16) 


(B— 17) 


It follows that the film-coefficient equation for mass transfer in the multi- 
component, chemically reacting boundary layer is the same as Eq. (B-14) with 
the Z* replaced by Z£ 


j k = p e u e C M<% < B - 18 > 

w e 

where is the mass transfer rate of chemical element k, independent of 

molecular configuration. 

Equation (B-18) is the desired form for a film coefficient approach to 
mass transfer in the chemically reacting, multicomponent boundary layer. It 
includes a driving potential that has been rationalized on a phenomenological 
basis. in the next section an approach for evaluating the mass- transfer co- 
efficient, P e u e C M , is described. 
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B.2.2 The Mass-Transfer Coefficient 

The mass- transfer "film-coefficient” is expressed as the product of the 

boundary-layer-edge mass velocity, p u , and a Stanton number for mass trans- 

e e 

fer, C w . It seems reasonable that the Stanton number for mass transfer in 
the multicomponent boundary layer is related to the Stanton number for heat 
transfer, C u , in the same manner as for a binary mixture 

ri 

C M = C H (Le) 3 / 3 (B-19) 

where the Lewis number is defined in terms of the effective diffusion coef- 
ficient defined earlier (Eq. (B-9) ) 

*# e f f 

Le = (B-20) 

with a the thermal diffusivity. Since the effective diffusion coefficient 

is a property of the mixture at a point, it will generally vary throughout 

the boundary layer and some boundary-layer reference state should be employed 

for its evaluation as is often done for the evaluation of all properties in 

the compressible boundary layer. For low-speed flows, the reference state 

could be taken as the arithmetic average between wall and boundary- layer 

edge conditions, whereas for high-speed flows, a reference state such as 

B-4 

that described by Eckert would seem approprxate. 
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APPENDIX C 


A REVIEW OF THE LARANJEIRA MODEL FOR MULT I COMPONENT 
THERMAL DIFFUSION COEFFICIENTS AS IT RELATES 
TO THE PRESENT MODEL 


In Ref. C-l, Laranjeira presents a theory for thermal diffusion factors 
based on free path concepts. On the basis of his physical model, he then 

postulates a first-order approximation to his theory which bears a striking 

T 

resemblance to the approximate model for D^ developed in the present re- 
port. A thorough evaluation of the rather complex dissertation of Laran- 
jeira has yet to be accomplished. However, at the present time it can be 
said that the Laranjeira model tends to substantiate the major assumptions 
which were made in the development of the present model, namely: 1) the 

proposed equation for D^ (Eq. (66)), and 2) the generalization of Eqs. 

(81) and (82) for a binary system to Eqs. (83) and (84) for a multicomponent 
system. His results suggest that a physically more sound correlation for 
thermal diffusion data might be = 1/(F^)^ rather than = 1/F^. How- 

ever, until further data correlation reveals that the former is indeed more 
accurate, the latter is currently favored since it leads to substantially 
simpler results. 

C-l 

The free path theory of Laranjeira is an extension of the earlier 

C-2 C-3 •• C-4 C- 5 

work of Frankel, Furry, Furth, and Whalley and Winter. It is 

based on two mean free paths: one, A being the mean distance over which 

the molecules of kind i should carry the attribute of the number density; 

and the other, being the mean distance to carry the attribute of the 

root-mean square of the thermal velocities. He postulates that these are 

related by 

= (1 + a ± ) A ± (C-l) 

where a^ would be expected to depend upon the force laws acting between 
colliding molecules. Specifically, he points out that a^ can be thought 
of as depending upon the "hardness" of molecular interactions. He investi- 
gated a rather extensive set of binary thermal diffusion data for the rare 
gases and found that a^ usually fell between 0.10 and 0.40. 

Although free path theory is mathematically inexact (see Ref. 6, p. 106) , 
the equation for diffusion velocities is shown for binary mixtures to be 
formally identical to that given by the exact Chapman-Enskog theory when 
expressed in terms of the coefficients for concentration diffusion, thermal 
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diffusion, and pressure diffusion. The inexactness enters from the fact 
that these coefficients are expressed in terms of mean free paths and thus 
are not identical with the exact form obtained by evaluation of the velocity 
distribution function. He does, however, use a highly refined form for the 
mean free path (nearly equivalent to Maxwell's mean free path corrected for 
persistence of velocities after collision, see Ref. 6, Chap. 5) . Larenjeira 
points out that this form has been satisfactory when comparing experiments 
with Chapman-Enskog 1 s theory. 

At this point, his equation offers little if any simplification to rigor- 
ous theory for multicomponent thermal diffusion. However, he then employs a 
first-order approximation to the l ^ 

= (C-2) 

where K is a constant for a gas mixture 


K 


1.16 TT 


n 






(C-3) 


and 


n 



(C-4) 


where n ^ is the number density of species j. Substituting into his more 
rigorous equations, he obtains the following result for binary mixtures. 


where 


[°ab]i 

1 - a ab A b~ K 

(C-5) 

2 x a A a + *b*b 

a T_ 

ab 

= (a a + ap/2 

(C-6) 

A « 

a 

\ k “a 

(C — 7 ) 


1 - “b A 

(C-8) 

A a " 

1 - a ab a 
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U) »< a - \ and a a - cr b 

(ii) Vl/7/1 — ► 00 and cr , /cr — ► 00 
oa Da 

At intermediate conditions, he shows that the effect of the approximation of 
^i suc ^ as yi e id a somewhat different concentration dependence. 

Laranjeira then extends his elementary theory to multicomponent mixtures, 
developing general expressions in terms of i k and a k (again using Eq. 

(C— 1) ) . Following the same procedure as for binary mixtures, he obtained 
the result 


55 



(C— 14) 



i -E 


k . ._ . .. El 


A. 

l 


a . 

l 


x jA 


a . 

— 1 


If the a^ are again considered as a universal constant, Eq. (C-14) becomes, 
with the aid of Eqs. (C-10) and (C-ll) 



(C-15) 


The [a ±j ] is related to , a multicomponent thermal diffusion ratio which 

he defined by 1 

V = I x i x k a ik (c " 16) 

1 k 

T 

The relation between his multicomponent and has not been studied 

at this time. It is significant to note, however, that the form of Eq. (C-15) 
is such that it does tend to corroborate the generalization from Eq. (82) to 
(84) which, although logical, was done without physical justification. 

Laranjeira generated thermal diffusion data for several isotopic ternary 
mixtures. He used this data in an attempt to substantiate his multicomponent 
theory with moderate success. To quote his general conclusions, "In spite of 
the elementary character of our treatment and the approximations which have 
been introduced in order to obtain simple equations, the most important fea- 
tures of thermal and pressure diffusion are satisfactorily explained and the 
agreement with experiments is reasonable." Further study of his correlations 
is required in order to assess the meaning of the word "reasonable." 
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